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Abstract 


Recently,  there  has  been  increasing  interest  from  the  research  community  in  multiple- 
input  multiple  output  (MIMO)  architectures  for  a  variety  of  radar  applications.  We 
are  interested  in  clarifying  the  angle  estimation  accuracy  of  coherent  or  co-located 
MIMO  radar  to  gain  an  accurate  understanding  of  its  potential  benefits.  The  theo¬ 
retical  two-way  antenna  radiation  patterns  are  derived  for  three  configurations  of  a 
linear  array:  a  phased  array  configuration,  a  MIMO-1  configuration  using  orthogonal 
waveform  transmission  on  all  elements,  and  a  MIMO-2  configuration  using  orthog¬ 
onal  waveform  transmission  on  the  two  end  elements.  Field  experimental  results 
obtained  with  an  8-element  X-band  linear  array  time-division  radar  and  trihedral 
targets  are  discussed.  Phased  array  radar  data  obtained  with  a  Mixed  Experiments 
and  Simulations  Approach  (MESA)  are  also  discussed.  It  is  found  that  the  experi¬ 
mental  mainbeam  patterns  match  the  theoretical  patterns.  MIMO-1  is  shown  to  have 
the  same  two-way  radiation  pattern  as  that  of  a  phased  array  radar  configuration. 
Compared  to  MIMO-1  and  phased  array,  MIMO-2  has  enhanced  angle  estimation 
accuracy,  lower  gain  and  higher  sidelobes,  while  only  requiring  two  orthogonal  wave¬ 
forms  on  transmit. 

Significance  for  defence  and  security 


Multiple-input  multiple  output  (MIMO)  radars  have  a  number  of  potential  benefits 
and  may  represent  a  future  evolution  of  current  phased  array  radars.  The  results  dis¬ 
cussed  in  this  report,  obtained  with  an  8-channel  X-band  linear  array  time-division 
MIMO  radar,  clarify  the  angle  estimation  accuracy  of  coherent  MIMO  radar.  This 
understanding  is  crucial  in  the  development  of  MIMO  radar  signal  processing  algo¬ 
rithms  and  applications. 
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Resume 


Recemment,  nous  avons  observe  un  interet  accru  de  la  part  de  la  communaute  de  la 
recherche  pour  les  architectures  a  entrees  multiples,  sorties  multiples  (MIMO)  pour 
une  variete  d’applications  radars.  Nous  voulons  clariher  l’exactitude  de  l’estimation 
de  bangle  des  radars  MIMO  coherents  ou  co-localises,  pour  obtenir  une  comprehen¬ 
sion  precise  de  leurs  benefices  potentiels.  Le  diagramme  de  rayonnement  de  l’antenne 
bidirectionnel  theorique  est  derive  pour  trois  configurations  d’un  reseau  lineaire  : 
une  configuration  avec  une  antenne  a  reseau  a  elements  en  phase,  une  configuration 
MIMO-1  qui  utilise  la  transmission  de  formes  d’onde  orthogonales  pour  tous  les  ele¬ 
ments,  et  une  configuration  MIMO-2  qui  utilise  la  transmission  de  formes  d’ondes 
orthogonales  aux  deux  elements  des  extremites.  Des  resultats  experimentaux  obtenus 
sur  1c  terrain  avec  un  radar  de  bande  X  a  8  canaux  et  a  repartition  dans  lc  temps  sont 
discutes.  Des  donnees  de  radar  a  reseau  a  elements  en  phase,  obtenus  avec  une  ap- 
proche  mixte  d ’experimentations  et  de  simulations  (MESA),  sont  egalement  discutes. 
On  y  trouve  que  les  faisceaux  principaux  des  diagrammes  de  rayonnement  experimen¬ 
taux  concordent  avec  ceux  des  diagrammes  theoriques.  On  y  montre  que  MIMO-1  a 
le  merne  diagramme  de  rayonnement  bidirectionel  que  la  configuration  radar  a  reseau 
a  elements  en  phase.  En  comparaison  de  MIMO-1  et  du  radar  a  reseau  a  elements  en 
phase,  MIMO-2  a  une  exactitude  accrue  de  l’estimation  de  Tangle,  un  gain  d’antenne 
plus  faible  et  des  lobes  secondaires  plus  importants,  par  contre  il  ne  necessite  que 
deux  formes  d’ondes  orthogonales  pour  la  transmission. 

Importance  pour  la  defense  et  la  securite 


Les  benefices  potentiels  des  radars  a  entrees  multiples,  sorties  multiples  (MIMO)  sont 
nombreux  et  ces  derniers  pourraient  representer  revolution  future  des  radars  a  reseau 
a  elements  en  phase  actuels.  Les  resultats  discutes  dans  ce  rapport,  obtenus  avec  un 
radar  MIMO  de  bande  X  a  8  canaux  et  a  repartition  dans  le  temps,  clarifient  l’exacti¬ 
tude  de  l’estimation  de  Tangle  des  radars  MIMO  coherents.  Cette  comprehension  est 
cruciale  pour  le  developpement  des  algorithmes  de  traitement  des  signaux  des  radars 
MIMO  et  lours  applications. 
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1  Introduction 


Recently,  there  has  been  increasing  interest  from  the  research  community  in  multiple- 
input  multiple  output  (MIMO)  architectures  for  a  variety  of  radar  applications  [1,  2, 
3].  The  MIMO  concept  was  first  developed  for  communications  systems,  for  which 
orthogonal  transmit-receive  channels  are  used  to  circumvent  channel  fading  [4],  In 
the  field  of  radar,  widely  separated  antenna  MIMO  radars  exploit  spatial  diversity, 
while  co-located  or  coherent  MIMO  radars  use  multiple  transmit  waveforms  to  obtain 
further  degrees  of  freedom  and  improve  radar  performance.  See  [5]  and  [6]  for  a 
comprehensive  literature  survey  of  MIMO  radar.  In  this  work,  attention  is  focused 
on  coherent  MIMO  radar,  which  requires  the  use  of  orthogonal  waveforms  on  transmit 
and  increased  signal  processing  on  receive,  and  therefore  has  increased  complexity  and 
flexibility  compared  to  phased  array  operation. 

Coherent  MIMO  radar  has  a  number  of  potential  benefits,  including  the  use  of  omni¬ 
directional  search  modes  [7],  enhanced  clutter  cancellation  [8],  and  increased  Doppler 
resolution  [9].  Because  of  the  increased  complexity  of  MIMO  radar,  it  is  important 
to  gain  an  accurate  understanding  of  its  potential  benefits.  In  particular,  it  has 
been  previously  stated  that  MIMO  radar  has  enhanced  angle  estimation  accuracy 
compared  to  phased  array  radar.  This  increased  accuracy  results  from  an  appar¬ 
ent  decrease  in  antenna  beamwidth  when  the  linear  array  is  operated  as  a  MIMO 
radar  [10,  11,  12,  13,  14], 

In  this  paper,  we  are  interested  in  clarifying  this  statement  about  the  enhanced  angle 
estimation  accuracy.  The  mathematical  formulation  to  support  our  analysis  is  pro¬ 
vided  in  Section  2.  We  then  provide  simple  but  yet  conclusive  experimental  results 
to  support  our  analysis.  The  experimental  methodology  is  described  in  Section  3  and 
the  results  are  presented  in  Section  4.  A  discussion  and  conclusions  are  provided  in 
Section  5. 
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2  Theoretical  antenna  patterns 


Two-way  antenna  radiation  patterns  are  derived  for  an  A-element  linear  array,  for 
uniform  phase  and  magnitude  illumination,  operating  in  three  distinct  configurations: 
Phased  Array,  MIMO-1,  and  MIMO-2.  In  the  Phased  Array  configuration,  a  single 
waveform  is  transmitted  coherently  across  the  array,  and  the  return  is  received  co¬ 
herently  on  all  elements.  For  the  MIMO-1  configuration,  each  element  transmits  a 
distinct  orthogonal  waveform,  and  the  return  is  received  coherently  on  all  elements. 
For  the  MIMO-2  configuration,  each  of  the  end  elements  transmits  a  distinct  orthog¬ 
onal  waveform,  and  the  return  is  received  coherently  on  the  first  IV  —  1  elements. 
See  Figure  1  for  an  illustration  of  the  configurations.  Each  element  in  the  array  has 
element  factor  E(9),  where  9  is  the  azimuth  angle.  A  linear  array  with  N  elements 
has  a  normalized  power  array  factor  An(9),  given  by 

sin2[Ar7r((i/A)  sinfl] 

N^}  N2  sin2[7r(d/A)sin0]’  W 


where  d  is  the  inter-element  spacing,  and  A  is  the  wavelength.  The  array  factor,  for 
narrow-band  and  far-held  approximations,  can  be  derived  by  summing  the  vector 
contributions  from  all  elements,  i.e. 


N- 1 

2 

A-n{9)  — 

-L  V  ^  A  27TC 

\  pJ  \  Tn 

VN  to 

1  N- 1 

= 

V''  j^nd  sine 

VNto 

(2) 

(3) 


where  rn  is  the  travel  time  between  a  point  target  and  the  nth  element  of  the  array. 


2.1  Phased  Array  configuration 

For  the  Phased  Array  configuration,  the  transmit  radiation  pattern  Gpa,tx($)  and 
the  receive  radiation  pattern  Gpa,rx(0)  are  equal,  so  that  [15] 

GW.Txt]#)  =  GpA,Rx(9)  =  E{9 )  An{9). 

The  Phased  Array  two-way  radiation  pattern  GpA,2-way(6l)  is  given  by 

G?PA,2-way(^)  =  Gta,Tx(0)  CpA,Rx(^)  =  E{9)2  A N(9)2 .  (4) 

2.2  MIMO-1  configuration 

The  signal  model  for  MIMO  operation,  for  which  M  orthogonal  waveforms  are  trans¬ 
mitted  with  uniform  amplitude  and  are  received  at  each  of  N  receive  elements,  for  a 
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Coherent  transmit  on  all  elements 


-Coherent  receive  on  all  elements - 

(a)  Phased  Array 


Orthogonal  transmit  on  all  elements 

i  i  i  1 

*  *  ‘  □  □ 


-Coherent  receive  on  all  elements - 

(h)  MIMO-1 


□ 


Orthogonal  transmit  on  end  elements  only 

□  □  *  *  *  □  □  □ 


Coherent  receive  on  first  N  —  1  elements 


(c)  MIMO-2 

Figure  1:  Configurations  for  a  linear  array  radar. 
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total  of  MN  matched  filters  and  MN  received  signals,  is  of  the  form: 


S 


M- 1 N- 1 


VMN  ,5o  S  e 


j¥(Tm+r„) 


1  M—l  /  1  iV-1 

_  ^  f  _L  V"  )  0i2¥Tm 


Vm  m=0  \y/N  n=0 


^  eJ  a  -™  j  e-*  a 


(5) 

(6) 


where  rm  and  rn  are  the  travel  times  between  a  point  target  and  the  mth  transmitter 
and  the  nth  receiver,  respectively.  This  equation  reduces  to 


S  =  An{6)Am{6) 


(7) 


when  considering  a  linear  array  with  the  narrowband  and  far-held  approximations. 

This  signal  model  applies  to  MIMO  beamforming  on  receive  and  is  equivalent  to  an 
effective  receive  array  TR(x)  that  is  the  spatial  convolution  of  a  transmit  array  T(x) 
and  a  receive  array  R(x)  [12].  Assuming 


T(x) 

M—l 

=  22  s(x  -  xm) 

777=0 

N-l 

(8) 

R(x) 

=  22^(X~Xn)’ 

n= 0 

(9) 

where  x  is  the  position  across  the  linear  array,  xm  are  the  locations  of  the  transmit 
elements  and  xn  are  the  locations  of  the  receive  elements,  then  the  spatial  convolution 
is  the  following: 


TR(x) 


Using 


then 


M- 1  N- 1 

T(x  —  s)R(s)ds  =  22  ~  xm  ~  s)  22  —  xn)ds 

m= 0  n= 0 


M—l  N—l  r 

22  22  -  (x  -  xm))5(s  -  xn)ds. 

771=0  71=0 


(10) 

(11) 


j  5(s- a)f(s)ds  =  /(a), 


(12) 


M—l  N- 1 

TR(x )  =  22  22  H(x  -  xm)  -  xn) 

771=0  n=0 
M—l  N-l 

=  22  22  HX~  (Xm  +  Xn)) 

771=0  71=0 


(13) 

(14) 
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which  is  also  a  linear  array  with  virtual  elements  located  at  x  —  xm  +  xn  for  0  <  m  < 
M  —  1  and  0  <  n  <  N  —  1. 

Using  this  signal  model  for  the  MIMO-1  configuration,  the  transmit  radiation  pattern 
Gmi,Tx(0)  is  given  by  the  element  pattern, 

GuiM0)  =  E(0)- 

On  receive,  the  effective  array  is  the  spatial  convolution  of  the  transmit  linear  array 
and  the  receive  linear  array.  Therefore,  the  MIMO-1  receive  array  factor  is  the  product 
of  an  N-element  linear  array  factor  by  itself,  that  is,  An{9)2.  The  MIMO-1  receive 
radiation  pattern  Gmi,Rx($)  is  then  specified  by 

GM1^(9)  =  E(9)AN(9)2. 

The  MIMO-1  two-way  radiation  pattern  GMi,2-way($)  is  given  by 

Cnv[l,2-way($)  =  GW,Tx($)  GW,Rx(#)  =  E{9)2  A n{9)2  .  (15) 

The  two-way  radiation  pattern  for  the  MIMO-1  configuration  is  equal  to  the  two-way 
radiation  pattern  for  the  Phased  Array  configuration. 

2.3  MIMO-2  configuration 

For  the  MIMO-2  configuration,  the  transmit  radiation  pattern  Gm2,Tx($)  is  given  by 
the  element  pattern, 

Gm2,Tx($)  =  E(6). 

On  receive,  the  effective  array  is  the  spatial  convolution  of  the  2-element  transmit 
linear  array  and  the  (N  —  l)-element  receive  linear  array,  i.e.  a  linear  array  with 
2 N  —  2  elements.  Therefore,  the  MIMO-2  receive  array  factor  is  that  of  a  linear 
array  with  2N  —  2  elements.  The  MIMO-2  receive  radiation  pattern  Gm2,rx(0)  is 
then  specified  by 

dnvi2,Rx($)  =  E(Q)  A2n- 2{9)- 

The  MIMO-2  two-way  radiation  pattern  G*m2, 2-way (9)  is  given  by 

Gm2, 2— way  (9)  =  Gm2,Tk  (9)GM2M0)  =  E(9)2A2N_2(9).  (16) 
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3  Experiment  methodology 


For  the  purpose  of  this  research,  a  time-division  8-channel  MIMO  radar  was  developed 
and  tested  on  an  open  field  using  trihedrals  as  targets.  The  targets  and  layout  of  the 
scenes  are  described  in  Section  3.1.  The  array  is  described  in  Section  3.2.  The  labo¬ 
ratory  equipment  assembled  and  used  to  serve  as  a  radar  is  described  in  Section  3.3. 
Data  was  also  collected  using  a  Vector  Network  Analyser  (VNA)  to  complement  an 
Agilent  ADS  simulation  to  mimic  a  phased-array  radar.  This  methodology,  named 
Mixed  Experiments  and  Simulations  Approach  (MESA),  is  described  in  Section  3.4. 

3.1  Scene 

The  experiment  was  conducted  in  an  open  field,  as  shown  in  Figure  2.  The  array  was 
mounted  on  a  tripod  approximately  2.35  m  above  the  ground.  One  trihedral  denoted 
PC,  of  69.6  cm  side,  was  positioned  at  35  m  range  from  the  center  of  the  array,  at 
an  angle  of  -4.8°  off  boresight  in  azimuth,  and  a  height  of  approximately  1.7  m.  This 
trihedral  acted  as  a  calibration  target  and  remained  in  the  scene  for  all  scenarios. 
One  or  two  trihedrals,  both  of  75  cm  side,  were  positioned  at  a  range  of  45  m  from 
the  array  center,  at  a  height  of  approximately  1.6  m,  except  for  one  scene  where  the 
trihedral  was  lowered  to  a  height  of  1.0  m.  These  trihedrals  were  used  as  targets,  and 
were  positioned  at  various  angles  off  boresight  depending  on  the  scenario.  They  are 
denoted  as  PI  to  P4  in  Table  1,  which  provides  an  overview  of  the  different  scenarios 
and  the  location  of  the  targets.  The  far-held  range  of  the  array  and  of  the  trihedrals 
is  26,  29  and  34  m  respectively,  according  to  Rff  =  2D2/X  at  9.0  GHz. 

Ground  truth  location  of  the  targets  and  array  for  this  experiment  was  obtained  using 
a  Leica  Total  Station  TS15  [16]. 

3.2  Array 

A  linear  array  was  assembled  out  of  10  Narda  640  horn  antennae  [17],  as  shown  in 
Figure  3.  The  8  central  elements  of  the  array  are  active  while  the  two  end  elements 
are  terminated  to  reduce  edge  effects.  The  antenna  elements  are  positioned  8.2  cm 
apart  and  are  operated  at  9.0  GHz.  At  this  frequency,  grating  lobes  are  expected 
every  24°. 

The  dimensions  of  the  Narda  640  horn  antennae  are  5.954  x  7.859  cm,  and  their 
beamwidths  are  estimated  to  be  30°  and  32°  at  9.0  GHz,  in  the  E  and  H  planes 
respectively.  For  this  experiment,  vertical  polarization  was  used. 
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Table  1:  List  of  runs  for  the  time-division  MIMO  experiment. 


F  ilename 

Description 

Target ( 

s) 

PC 

PI 

P2 

P3 

P4 

-4.8° 

-0.4° 

+1.6° 

+4.6° 

+14.8° 

35  m 

45  m 

45  m 

45  m 

45  m 

scene 1 

Empty  scene 

X 

scene5 

One  target  at  -0.4° 

X 

X 

scene7 

One  target  at  +1.6° 

X 

X 

scene3 

One  target  at  +4.6° 

X 

X 

scene2 

One  target  at  +4.6°,  lower 
height 

X 

X 

scene9 

One  target  at  +14.8° 

X 

X 

scene6 

Two  targets  2.1°  apart 

X 

X 

X 

scene4 

Two  targets  5.0°  apart 

X 

X 

X 

scene8 

Two  targets  13.2°  apart 

X 

X 

X 

Figure  2:  Photograph  of  the  experimental  setup,  showing  the  array,  the  smaller 
trihedral  used  as  a  calibration  target,  and  the  larger  trihedral  used  as  the  target. 
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Figure  3:  Photograph  of  the  linear  array  assembled  for  the  time-division  MIMO 
experiment. 
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3.3  Time-division  8-channel  MIMO  radar 


3.3.1  Measurement  setup 

A  time-division  8-channel  MIMO  radar  was  assembled  using  laboratory  equipment  as 
shown  in  Figure  4.  A  list  of  components  used  for  the  test  is  provided  in  Table  2.  An 
Arbitrary  Waveform  Generator,  (AWG1),  produces  differential  I/Q  signals  which  are 
routed  to  the  wideband  input  ports  of  a  Vector  Signal  Generator,  VSG1.  The  VSG 
upconverts  the  base-band  I/Q  signal  to  X-band.  This  results  in  an  Linear  Frequency 
Modulated  (LFM)  signal  centered  at  9.0  GHz,  with  a  bandwidth  of  150  MHz  and 
a  pulse  width  of  100  /zs.  The  chirp  rate  is  +1.5  MHz//zs.  The  Tx  signal  is  routed 
sequentially  to  each  of  the  8  Tx/Rx  antennas  via  a  SP10T  microwave  switch,  SW1. 
The  received  signal  is  routed  via  SP2T  switches,  SW2-SW5,  to  the  4  channels  of  a 
digital  oscilloscope,  OSC1.  The  measurement  is  fully  automated  and  uses  Agilent’s 
Signal  Studio  for  Pulse  Building  (SSPB),  N7620B,  to  generate  the  transmit  wave¬ 
form,  Agilent’s  Vector  Signal  Analysis  (VSA)  software  to  demodulate,  decimate,  and 
perform  data  transfer  of  the  digitized  time  domain  data.  Control  of  the  SSPB  and 
VSA  software  along  with  switch  control  and  data  storage  is  accomplished  via  in-house 
code  written  in  Agilent’s  VeePro  9.3  software.  The  VeePro  software  interfaces  with 
the  SSPB  and  VSA  software  packages  via  .COM  and  .NET  APIs  respectively. 

Prior  to  measurement,  a  calibration  is  performed  in  order  to  flatten  the  frequency 
response  of  the  Tx  signal.  The  spectrum  analyzer,  SAl,  is  connected  to  the  ends  of 
each  of  the  RF  cables  (CAB1-8)  in  turn.  The  AWG  and  VSG  produce  a  GW  tone 
which  is  swept  across  a  frequency  band  centered  at  9.0  GHz  with  a  bandwidth  of 
250  MHz.  The  SSPB  software  then  uses  this  information  to  predistort  the  base  band 
I/Q  signals  such  that  a  flat  frequency  response  is  obtained  at  the  antenna  input. 
The  calibrated  signal  for  each  of  the  8  paths  of  the  Tx  signal  are  stored  and  applied 
during  the  measurement.  During  the  calibration,  the  power  at  the  antenna  input  is 
measured  and  adjusted  such  that  it  is  equal  for  all  8  channels. 

The  measurement  procedure  involves  transmitting  a  signal  on  one  transmit  channel, 
then  receiving  it  on  the  4  odd  numbered  receive  channels  simultaneously  (i.e  receive 
channels  Rxl,  Rx3,  Rx5,  and  Rx7).  The  transmit  switch,  SW1,  is  then  switched  to 
the  next  transmit  channel  and  the  signal  is  received  on  the  same  4  odd  numbered 
receive  channels.  This  procedure  is  repeated  until  the  Tx  signal  has  been  transmitted 
through  all  8  Tx  channels.  The  receive  switches,  SW2-SW5,  are  then  toggled  such 
that  the  even  numbered  receive  channels,  Rx2,  Rx4,  Rx6,  and  Rx8  are  now  connected 
to  OSC1  and  the  odd  numbered  Rx  channels  are  terminated  into  50  O.  The  Tx  signal 
is  now  transmitted  in  sequence  through  each  Tx  channel  (via  SW1)  and  received  on 
the  4  even  numbered  Rx  channels  simultaneously. 

When  the  measurement  procedure  is  finished,  data  for  a  complete  8x8  time  division 
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Table  2:  List  of  components  for  the  time-clivision  MIMO  experiment. 


Designator 

Part 

Manufacturer 

Description 

ADI-8 

— 

— 

N(f)  to  SMA(m)  adapter 

AD9-12 

54855-67604 

Agilent 

Precision  BNC  to  3.5  mm  (f) 
adaptors 

ANTI¬ 
ANT  10 

640 

Narda 

8-12.4  GHz  Horn,  No.  5 

AWG1 

M8190A 

Agilent 

Arbitrary  Waveform  Generator 

12  Gsa/s 

VSG1 

E8267D 

Agilent 

PSG  Vector  Signal  Generator, 

20  GHz 

SA1 

E4440 

Agilent 

Performance  Spectrum  Analyzer 
(PSA) 

OSC1 

DSO91204A 

Agilent 

Digital  Oscilloscope,  BW  = 

12  GHz 

CAB  1-8 

SF104PEA/16N/11SMA 

Huber  &  Suhner 

2.5  m  rf  cables 

CAB9-16 

MiniBend-16 

Astrolab 

16"  rf  cables 

CAB17-20 

MiniBend-12 

Astrolab 

12"  rf  cables 

CAB21-24 

MiniBend-16 

Astrolab 

16"  rf  cables 

CAB25, 

CAB27, 

CAB28 

EF400 

Huber  &  Suhner 

1  m  BNC  cable 

CAB26 

— 

Huber  &  Suhner 

20"  BNC  cable 

CAB29 

— 

— 

GPIB  cable 

CAB30 

— 

— 

GPIB  cable 

CAB31 

EF400 

Huber  &  Suhner 

2  m  BNC  cable 

CIR1-CIR8 

ESC0950 

Star  Microwave 

8-11  GHz  Microwave  Circulator 

ATI 

MCL  BW 

Mini-Circuits 

dc-18  GHz  3  dB  attenuator 

AT2 

FP-50 

Texscan 

3  dB  BNC  attenuator 

Tl-2 

MCL  ANNE-50+ 

Mini-Circuits 

dc-20  GHz  50  fi  termination, 
SMA(m) 

DC1,  DC2 

E3648A 

Agilent 

100  W  dual  output  DC  power 
Supply 

DC3 

E3620A 

Agilent 

50  W  dual  output  DC  power 
Supply 

DC4 

E3649A 

Agilent 

100  W  dual  output  DC  power 
Supply 

SP1 

T1000 

Anzac 

10-1,000  MHz  power  splitter 

SW1 

MSN-  10DT-05-DEC-MP 

AMC 

SP10T  absorptive  switch, 

0.5  GHz-18  GHz 

SW2-5 

8762B  TTL  T24 

Agilent 

SP2D  coaxial  switch,  dc-18  GHz 

SW6 

DSS-8+ 

D-Link 

8  10/100  port  LAN  switch 

GPIB1 

82357B 

Agilent 

USB  to  GPIB  interface 

PCIE1 

M9045B 

Agilent 

PCIe  laptop  card  adapter  Gen 

1x4 

PCI 

HP  Elitebook  8560P 

Hewlett  Packard 

Laptop  running  Vee  Pro  9.3, 

VSA  16.0,  Signal  Studio  for 

Pulse  Building 
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MIMO  system  is  obtained.  The  data  is  stored  in  Matlab  compatible  .mat  files. 

3.3.2  Signal  processing  algorithms 

The  signal  processing  algorithms  consist  of  pulse  compression,  array  calibration,  and 
MIMO  beamforming. 

De-chirping  of  the  64  received  signals  is  performed  by  complex  multiplication  with 
the  complex  conjugate  of  the  transmitted  waveform.  A  fast  Fourier  transform  then 
provides  a  pulse  compressed  range  profile  for  each  of  the  64  channels. 

Array  calibration  is  performed  by  using  the  64  range  profiles  and  the  known  positions 
of  the  antenna  array  elements  and  of  the  calibration  target  located  at  35  m  range. 
For  each  transmit-receive  element  pair,  the  delay  and  phase  due  to  propagation  from 
the  transmit  element  to  the  calibration  target  and  back  to  the  receive  element  are 
determined.  The  corresponding  range  profile  is  then  shifted  in  range.  This  is  a 
coarse  correction  for  cable  lengths  and  various  delays  introduced  by  the  equipment. 
A  finer  correction  is  then  applied  as  a  phase  correction.  This  process  is  repeated 
for  all  combinations  of  transmit  and  receive  channels,  producing  64  calibrated  range 
profiles. 

MIMO  beamforming  is  performed  using  the  following  delay-and-sum  method.  The 
two-way  distance  and  phase  due  to  propagation  from  transmit  antenna  A  to  a  given 
point  in  space  and  back  to  receive  element  B,  is  determined.  The  two-way  distances 
and  phases  for  all  combinations  of  transmit  and  receive  elements  are  determined. 
The  two-way  distances  are  used  to  pick  the  signals  as  a  linear  interpolation  from 
the  appropriate  range  bins  in  the  calibrated  range  profiles.  The  phases  are  used  to 
compensate  the  received  signals.  For  beamforming  using  a  MIMO-1  configuration, 
the  resulting  64  phase  compensated  complex  signals  corresponding  to  the  point  in 
space  are  summed.  For  beamforming  using  any  other  MIMO  configuration  such  as 
MIMO-2,  the  appropriate  subset  of  the  64  phase  compensated  complex  signals  are 
summed.  This  process  is  repeated  for  all  points  in  the  space  of  interest.  Typically, 
a  grid  in  space  is  defined  at  the  range  of  the  targets  (around  45  meters)  and  for 
angles  varying  between  -40  and  +40°.  In  Section  4,  the  results  are  plotted  as  an 
azimuth  profile  at  the  range  of  maximum  intensity  as  a  function  of  azimuth  angle. 
The  azimuth  profile  is  displayed  as  relative  power  levels  or  as  normalized  power,  in 
units  of  decibels.  This  azimuth  profile  indicates  the  beamforming  capability  of  the 
considered  array. 

Given  that  the  MIMO  beamforming  algorithm  for  this  report  is  a  simple  sum  of  the 
appropriate  signals  and  no  scaling  is  applied,  we  expect  different  power  levels  between 
the  MIMO  configurations.  Explicitly,  the  power  level  in  decibels  scales  as 

101og10(iVTxiVRx)2,  (17) 
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where  NTx  and  NRx  are  the  number  of  transmit  and  receive  elements,  respectively. 
For  MIMO-1,  NTx  =  NRx  =  8  and  for  MIMO-2,  Ntx  =  2  while  NRx  =  7,  leading  to 
a  13  dB  difference. 


12 


DRDC-RDDC-201 5-R051 


ANTI  ANT2  ANT3  ANT4  ANT5  ANT6  ANT  7  ANT  8  ANT9  ANT10 


experiment. 
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3.4  Mixed  Experiments  and  Simulations  Approach  (MESA) 


The  purpose  of  MESA  is  to  determine  the  frequency  response,  and  hence  the  impulse 
response  of  the  system  composed  of  the  8  antennae  and  scene  of  observation  including 
the  targets.  By  determining  this  information,  the  data  can  be  used  with  a  system 
simulator  to  study  the  effects  of  changing  system  parameters  such  as  Tx  waveform 
characteristics  and  adding  impairments  such  as  nonlinear  distortion  and  noise  to  the 
system.  In  this  study,  the  8-port  S-parameters  of  the  antenna  array  and  observation 
scene  were  measured  and  the  data  was  imported  into  Agilent  ADS  electronic  design 
automation  (EDA)  software  in  order  to  simulate  an  ideal  phased-array  radar.  The 
measurement  and  simulation  procedure  is  described  in  the  following  sections. 

3.4.1  Measurement  Setup 

The  test  setup  used  to  measure  the  S-parameters  is  shown  in  Figure  5  and  the  equip¬ 
ment  list  is  described  in  Table  3.  The  VNA  is  connected  to  the  12  port  test  set  such 
that  the  full  8-port  S-parameters  of  the  antenna-array  and  scene  can  be  measured.  A 
two-port  ECAL  module  along  with  an  unknown  THRU  standard  is  used  to  perform 
the  8-port  calibration  [18].  The  unknown  THRU  consists  of  a  high  quality  phase 
stable  microwave  cable  and  two  SMA  adapters.  One  port  of  the  E-CAL  module  is 
connected  to  each  port  of  the  test  set  and  then  the  unknown  THRU  is  successively 
connected  between  two  ports  of  the  test  set.  The  VNA  software  minimizes  the  num¬ 
ber  of  connections  necessary  between  the  ports  in  order  to  extract  all  the  error  terms 
required  for  the  calibration. 

The  S-parameters  were  measured  for  run  scene5,  consisting  of  PI  and  PC  and  run 
scene9,  consisting  of  P4  and  PC. 

3.4.2  VNA  Setup 

In  order  to  obtain  a  good  impulse  response,  the  VNA  is  setup  to  sweep  over  a  wide 
band.  The  lower  frequency  is  limited  by  the  VNA  at  10  MHz  and  the  higher  frequency 
is  set  to  18  GHz,  which  is  the  highest  frequency  for  which  the  components  in  the  setup 
are  rated. 

The  frequency  step  size  must  be  small  enough  such  that  the  frequency  response  is  not 
undersampled.  The  change  in  phase  A0  of  the  S-parameter  per  change  in  frequency 
A /  must  be  less  than  180°  based  on  the  two-way  time  delay  to  the  furthest  target  of 
interest,  i.e. 


A  0  = 
or  A  /  < 


(18) 

(19) 
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VNA1 


Observation  Scene 


Figure  5:  Setup  used  to  measure  multi-port  S-parameters 
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Table  3:  Equipment  used  for  S-parameter  measurement. 


Designator 

Part 

Manufacturer 

Description 

ADP1-8 

- 

~ 

N-female  to  SMA-male 
adapter 

ADP9-10 

Southwest 

Microwave 

SMA-female  to 
SMA-female  adapter 

ANTI-10 

640 

Narda 

8.2  -  12.4  GHz  standard 
gain  horn 

CALI 

N4691-60006 

Agilent 

300  kHz  -  26.5  GHz 

E-cal  module 

CBL1-9 

Sucoflex 

104PEA 

Huber  &  Suhner 

2.5  m  RF  cables  N-male, 
SMA  male 

CBL10 

Accu-Test  150 

Teledyne  Storm 

18"  RF  cable 

TRM1-2 

ANNE-50X+ 

Mini- Circuits 

50  SMA-male 
termination 

VNAl 

N5242A 

Agilent 

4  port  VNA 

VNA2 

U3042A  E12 

Agilent 

12  port  test  set 

P4 

DREO 

— 

Trihedral  corner 
reflector,  l  =  75  cm 

PI 

DREV 

— 

Trihedral  corner 
reflector,  l  =  75  cm 

PC 

Athena 

Trihedral  corner 
reflector,  /  =  69.6  cm 
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where  Rmax  is  the  maximum  range  of  interest  and  c  is  the  speed  of  light.  Assuming 
Rmax  =  50  m,  then  it  is  found  that  the  frequency  step  A /  must  be  less  than  1.5  MHz. 
However,  for  the  software  to  calculate  a  good  impulse  response  it  is  suggested  that 
there  should  be  at  least  6  samples  per  360°  phase  change  [19].  As  a  result,  (19)  is 
modified  to: 


A/< 


c 


10/L 


(20) 


and  A /  <  600  kHz. 


Selection  of  the  IF  bandwidth,  BWif  involves  a  tradeoff  between  measurement  speed 
and  noise  floor.  Based  on  the  PNA-X  specifications  [20],  the  typical  noise  floor  of  the 
test  ports  at  X-band  is  -117  dBm  in  a  10  Hz  bandwidth  which  translates  to  a  noise 
density,  Nd  of  -127  dBm/Hz.  The  noise  floor  of  the  measurement  can  be  calculated 
from: 

N floor  =  NdBWlF .  (21) 


An  estimate  of  the  received  signal  strength  can  be  made  using  the  radar  range  equa¬ 
tion  and  the  ideal  radar  cross  section  for  a  trihedral  corner  reflector.  For  a  trihedral 
corner  reflector  with  sidelength  l,  at  wavelength  A  the  radar  cross  section,  crtgt  is  given 
by  [21]: 


AttB 
<T"" =  3W 


(22) 


The  received  signal-to-noise  ratio  is  given  by: 


SNR  = 


PrGlnt  A2 
(47 TfR±Nfloor 


atgt 


(23) 


where  Gant  is  the  gain  of  the  antenna,  R  is  the  range,  and  Pt  is  the  transmit  power. 


Substituting  (21)  and  (22)  into  (23)  results  in 

Gant 
An 

Assuming  the  following: 


SNR  = 


3NdBWjF 


(24) 


Pt  is  the  transmit  power,  8  dBm.  This  is  the  maximum  power  that  could 
be  transmitted  without  overloading  the  instrument. 

Nd  is  the  noise  density  of  the  VNA  input  ports,  -127  dBm/Hz  typical  at 
X-band. 

BWif  is  the  bandwidth  of  the  VNA  IF  receiver.  This  parameter  is  adjustable 
and  was  set  to  1  kHz. 

R  is  the  range  to  the  farthest  target,  45  m. 

I  is  the  side  length  of  the  corner  reflector,  75  cm. 

Gant  is  the  gain  of  the  antenna,  15.4  dBi. 
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The  SNR  is  approximately  38  dB  which  was  more  than  adequate  for  this  measure¬ 
ment.  A  summary  of  all  the  VNA  settings  for  the  measurement  is  given  in  Table  4. 

Table  4:  VNA  Parameters  used  for  measurement  and  calibration. 


Parameter 

Value 

Units 

Start  Frequency 

10 

MHz 

Stop  Frequency 

18.01 

GHz 

Step  Frequency 

562.5 

kHz 

No.  of  points 

32001 

— 

IF  Bandwidth 

1 

kHz 

Transmit  Power 

8 

dBm 
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3.4.3  Simulation 


The  measured  8-port  S-parameters  for  the  observation  scene  were  imported  into  Ag¬ 
ilent  ADS  in  order  to  simulate  an  ideal  phased-array  radar.  The  first  step  was  to 
calibrate  the  array  to  account  for  differences  between  the  antenna  ports  due  to  cable 
lengths,  adaptors,  etc.  This  was  done  by  performing  a  transient  simulation.  Each 
port  of  the  array  was  independently  stimulated  with  an  linear  frequency-modulated 
(LFM)  signal  at  a  center  frequency  of  9  GHz  and  bandwidth  of  2  GHz.  The  monos¬ 
tatic  returns  from  each  port  were  multiplied  by  the  replica  of  the  transmit  signal  and 
the  resulting  signal  was  transformed  to  the  frequency  domain.  The  returns  from  each 
target  appear  at  frequencies  proportional  to  their  range  and  the  chirp  rate.  A  result 
is  shown  in  Figure  6(a)  for  an  uncalibrated  array.  The  appropriate  delays  were  added 
to  each  port  of  the  array  until  all  the  peaks  were  aligned  and  the  phase  differences 
between  the  ports  were  minimized.  This  result  is  shown  in  Figure  6(b). 

The  phased-array  simulations  were  carried  out  using  an  LFM  source  at  a  center 
frequency  of  9  GHz  with  a  bandwidth  of  150  MHz,  the  same  parameters  as  were  used 
for  the  MIMO  measurements,  as  described  in  Section  3.3.  The  array  was  steered  over 
the  observation  scene  at  a  determined  angle  and  the  returns  were  pulse  compressed 
and  then  power  combined  to  obtained  the  final  output.  For  example  using  the  data 
from  run  scene9,  when  the  array  is  steered  to  14.8°  such  that  it  is  pointing  directly 
at  P4,  we  obtain  the  output  power  as  a  function  of  range  shown  in  Figure  7.  The 
returns  from  P4,  PC  and  also  those  due  to  antenna  return  loss  are  clearly  visible. 

For  each  data  set,  the  array  was  steered  to  angles  ±15°  around  the  target  of  interest, 
and  the  resulting  target  output  power  was  determined  as  a  function  of  angle. 

Alternatively,  the  measured  8-port  S-parameters  were  processed  as  a  MIMO  array. 
The  same  LFM  source  was  used  in  the  Agilent  ADS  simulations  but  the  array  was 
not  steered  and  the  channels  were  not  power  combined.  The  pulse  compressed  range 
profiles  were  processed  as  described  in  Section  3.3.2. 
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Figure  6:  Returns  from  calibration  target  PC.  The  8  colours  correspond  to  the  8 
receive  channels. 
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Figure  7:  Received  signal  from  phased  array  simulation  for  scene  9  at  azimuth  angle 
=  14.8?  The  antenna  array  is  electrically  steered  toward  target  P4. 
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4  Experimental  results 


In  this  Section  are  presented  the  experimental  results  obtained  with  the  MIMO  ar¬ 
ray.  The  results  obtained  with  the  time-division  8-channel  radar,  processed  with  a 
MIMO-1  configuration,  are  presented  in  Section  4.1.  The  MESA  data,  collected  us¬ 
ing  the  VNA  and  processed  in  simulations  using  the  ADS  software,  are  presented  in 
Section  4.2  to  demonstrate  that  the  MIMO-1  and  phased  array  factors  are  identical. 
Finally,  a  comparison  between  the  MIMO-1  and  MIMO-2  configurations  is  provided 
in  Section  4.3  using  the  time-division  8-channel  radar  data. 


4.1  MIMO-1  array  (8  transmit  elements  and  8  receive  elements) 

In  this  section  are  grouped  the  results  obtained  by  processing  the  time-division  radar 
data  using  a  MIMO-1  configuration  with  8  transmit  and  8  receive  elements. 


4.1 .1  One  target  scenes 


In  Figures  8  to  12  are  shown  the  experimental  data  beamformed  with  a  MIMO-1 
configuration  for  the  scenes  with  one  target.  The  MIMO  beamforming  algorithms 
were  described  in  Section  3.3.2.  The  theoretical  array  factors  for  the  2-way  radiation 
pattern  are  shown  for  comparison,  ignoring  the  element  factor.  The  theoretical  array 
factors  are  scaled  in  amplitude  to  match  the  experimental  data  and  steered  in  azimuth 
toward  the  target  location  at  90,  according  to: 


sin2  [AOr(<i/A)  (sin  9  —  sin0o)] 
N2  sin2[7r(cf/A)(sind  —  sin6*0)]  ’ 


(25) 


For  each  of  these  graphs,  we  see  that  the  main  beam  matches  the  ground  truth  of 
the  target  location.  We  observe  the  grating  lobes,  as  expected.  The  experimental 
data  matches  the  theoretical  antenna  pattern  especially  for  the  main  lobe  and  the 
first  sidelobe.  Scene  2  (Figure  11)  matches  particularly  well  the  theoretical  antenna 
pattern,  even  for  the  sidelobes. 

The  experimental  data  for  the  same  five  scenes  are  overlaid  on  the  graph  of  Figure  13, 
to  illustrate  the  difference  in  power  levels.  The  Figure  shows  the  azimuth  profiles 
after  beamforming,  but  the  range  compressed  responses  exhibit  the  same  power  level 
variations  from  scene  to  scene  (graph  not  shown).  These  power  level  variations  cannot 
be  explained  by  the  element  factor  of  the  antenna  element,  which  has  a  beamwidth 
of  approximately  32°  in  the  azimuth  direction  [17].  The  element  factor  accounts  for 
only  a  6  dB  decrease  between  the  target  at  boresight  and  the  one  located  at  +14.8°. 
Even  scenes  2  and  3,  with  the  target  located  at  the  same  position  in  azimuth,  show 
a  7  dB  difference  in  the  main  beam  power  level.  In  Annex  A,  it  is  demonstrated 
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Figure  8:  Azimuth  profile  for  scene  5  with  one  target  located  at  -0.4°.  The  theoret¬ 
ical  array  factor,  scaled  and  steered  in  azimuth,  is  shown  for  comparison. 


that  ground  multipath  can  most  likely  explain  the  different  power  levels,  especially 
for  scenes  2  and  3  with  the  target  at  the  same  azimuth  location  but  different  heights 
above  the  ground. 
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Figure  9:  Azimuth  profile  for  scene  7  with  one  target  located  at  +1.6°.  The  theo¬ 
retical  array  factor,  scaled  and  steered  in  azimuth,  is  shown  for  comparison. 


Figure  10:  Azimuth  profile  for  scene  3  with  one  target  located  at  +4.6°.  The 
theoretical  array  factor,  scaled  and  steered  in  azimuth,  is  shown  for  comparison. 


24 


DRDC-RDDC-201 5-R051 


Figure  11:  Azimuth  profile  for  scene  2  with  one  target  located  at  +4.6°,  lower  height. 
The  theoretical  array  factor,  scaled  and  steered  in  azimuth,  is  shown  for  comparison. 


Figure  12:  Azimuth  profile  for  scene  9  with  one  target  located  at  +14.8°.  The 
theoretical  array  factor,  scaled  and  steered  in  azimuth,  is  shown  for  comparison. 
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Figure  13:  Azimuth  profile  for  scenes  5,  7 ,  3,  2,  and  9  with  one  target  at  various 
locations  in  azimuth. 
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4.1 .2  Two  target  scenes 


In  Figures  14  to  16  are  shown  the  time-division  radar  data  for  the  scenes  with  two 
targets.  For  each  graph,  the  experimental  data  of  the  corresponding  scenes  with  one 
target  only  are  displayed  for  comparison. 


Figure  14:  Azimuth  profile  for  scene  6  with  two  targets  located  at  -0.4°  and  + 1.6 °. 
The  azimuth  profiles  for  scene  5  (one  target  at  -0.4°)  and  scene  7  ( one  target  at  +1.6°) 
are  displayed  for  comparison. 


For  scene  6  with  two  targets  located  2.1°  apart  (Figure  14),  the  resolution  of  the 
MIMO-1  array  is  not  sufficient  to  discriminate  between  the  two  targets. 

For  scene  4  with  two  targets  located  5°  apart  (Figure  15),  the  two  targets  are  resolved, 
but  we  notice  that  the  target  at  +4.6°  (scene  4,  red  curve)  appears  with  a  relative 
power  1.3  dB  lower  than  it  appears  when  it  is  the  only  target  in  the  scene  (scene  2, 
green  curve).  On  the  other  hand,  the  target  at  -0.4°  (scene  4,  red  curve)  appears  with 
a  relative  power  1.5  dB  higher  than  it  appears  when  it  is  the  only  target  in  the  scene 
(scene  5,  blue  curve).  The  target  at  +4.6°  was  not  moved  between  the  scene3  and 
scene4  data  acquisitions,  and  the  target  at  -0.4°  was  not  moved  between  the  scene4 
and  scene5  data  acquisitions.  It  is  hypothesized  that  there  is  interaction  between  the 
two  trihedrals  on  the  field. 

For  scene  8  with  two  targets  located  13.2°  apart  (Figure  16),  the  two  targets  are 
resolved.  The  two  trihedrals  are  far  enough  to  not  interact,  and  the  experimental 
data  for  two  targets  coincide  with  the  experimental  data  for  each  target  individually, 
as  expected. 
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Figure  15:  Azimuth  profile  for  scene  4  with  two  targets  located  at  -0.4°  and  +4.6°. 
The  azimuth  profiles  for  scene  5  (one  target  at  -0.4°)  and  scene  3  (one  target  at  +4.6°) 
are  displayed  for  comparison. 


angle  (degrees) 

Figure  16:  Azimuth  profile  for  scene  8  with  two  targets  located  at  +1.6°  and  +14.8°. 
The  azimuth  profiles  for  scene  7  (one  target  at  +1.6°)  and  scene  9  (one  target  at 
+14.8°)  are  displayed  for  comparison. 
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4.2  Phased  array  results  obtained  with  MESA  (8  transmit 
elements  and  8  receive  elements) 

In  this  section,  it  is  demonstrated  experimentally  that  the  phased  array  2- way  antenna 
pattern  and  the  MIMO-1  configuration  have  identical  array  factors.  In  Figures  17 
and  18  are  shown  the  results  for  scenes  9  and  5,  respectively.  Each  graph  shows  four 
curves.  The  theoretical  two-way  array  factor  is  shown  as  a  solid  line.  The  time- 
division  radar  data  processed  with  a  MIMO-1  configuration  are  shown  as  dots.  The 
MESA  data  processed  as  a  phased  array  simulation,  as  described  in  Section  3.4,  are 
shown  as  circles.  Finally,  the  MESA  data  processed  as  a  MIMO-1  configuration  are 
shown  as  triangles.  All  curves  are  shown  on  a  normalized  scale  since  the  time-division 
radar  data  and  the  MESA  data  have  different  received  power  levels  due  the  different 
setup  and  cables. 


angle  (degrees) 


Figure  17:  Azimuth  profile  for  scene  5  with  one  target  located  at  -0.4°,  for  the 
time-division  MIMO-1  radar  data,  the  MESA  data  processed  as  phased  array  and  as 
MIMO-1.  The  theoretical  array  factor,  steered  in  azimuth,  is  shown  for  comparison. 


The  four  curves  are  in  agreement  in  terms  of  the  main  beam  location  and  width.  We 
see  a  discrepancy  in  the  sidelobes.  The  MESA  data,  either  processed  as  MIMO-1  or 
phased  array,  show  lower  siclelobe  levels  than  the  time-division  radar  data.  The  fact 
that  the  MESA  data  shows  similar  sidelobe  levels  for  both  MIMO-1  and  phased  array 
simulations  indicates  that  the  higher  sidelobes  of  the  time-division  data  are  not  due 
to  the  MIMO  processor.  Rather,  the  time-division  data  were  acquired  as  8-bit  data 
and  suffer  from  lower  dynamic  range.  Furthermore,  they  were  acquired  with  multiple 
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Figure  18:  Azimuth  profile  for  scene  9  with  one  target  located  at  +14.8°,  for  the 
time-division  MIMO-1  radar  data,  the  MESA  data  processed  as  phased  array  and  as 
MIMO-1.  The  theoretical  array  factor,  steered  in  azimuth,  is  shown  for  comparison. 


triggers  which  can  cause  jitter  issues.  We  believe  that  these  effects  are  small  enough 
to  not  affect  the  study  of  the  MIMO  angle  accuracy,  but  are  large  enough  to  affect 
the  sidelobe  structure  levels. 
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4.3  Comparison  of  MIMO-1  and  MIMO-2  arrays 


In  this  section  are  compared  the  azimuth  profiles  of  the  time-division  radar  data 
processed  using  both  MIMO-1  and  MIMO-2  arrays. 

4.3.1  One  target  scene 

In  Figure  19  is  shown  the  time-division  radar  data  for  scene  2  with  one  target  at 
+4.6°  and  lower  height,  processed  for  both  MIMO-1  and  MIMO-2  arrays.  On  the  top 
graph,  the  azimuth  profiles  are  shown  as  normalized  values.  We  can  clearly  see  the 
improved  resolution  of  the  MIMO-2  array  over  the  MIMO-1  array,  at  the  expense  of 
larger  sidelobes  (13  dB  versus  26  dB  main-to-first  sidelobe  ratio  for  the  MIMO-2  and 
MIMO-1  arrays,  respectively).  On  the  bottom  graph,  the  azimuth  profiles  are  shown 
without  normalization  and  we  see  the  gain  loss  of  approximately  13  dB  for  MIMO-2 
array,  as  expected  since  only  two  transmit  elements  are  used  instead  of  all  8  elements 
for  the  MIMO-1  array. 

4.3.2  Two  target  scenes 

In  Figures  20  to  22  are  shown  the  azimuth  profiles  for  the  scenes  with  two  targets, 
processed  for  both  MIMO-1  and  MIMO-2  arrays.  We  see  the  effect  of  slightly  in¬ 
creased  resolution  for  all  cases,  at  the  expense  of  larger  sidelobes  for  MIMO-2.  Given 
the  issues  of  ground  multipath  and  interaction  between  targets  discussed  above,  the 
resulting  signatures  for  the  two  targets  can  result  in  power  values  differing  by  as  much 
as  13  dB.  Using  a  MIMO-2  array,  the  weaker  target  now  appears  at  approximately  the 
same  power  level  as  the  sidelobes  of  the  stronger  target.  As  a  result,  the  signature  of 
the  weaker  target  appears  as  the  constructive  or  destructive  sum  of  the  target  signal 
itself  and  the  stronger  target  sidelobes.  In  particular,  the  weaker  target  of  scene  8 
appears  10  dB  lower  than  the  stronger  target  for  MIMO-1,  while  it  appears  14  dB 
lower  for  MIMO-2  (Figure  22). 


DRDC-RDDC-201 5-R051 


31 


received  power  (dBm)  normalized  power  (dB) 


(b)  Received  power  values 

Figure  19:  Azimuth  profiles  for  scene  2  with  one  target  located  at  +4.6°,  lower 
height,  processed  for  both  MIMO-1  and  M1MO-2  arrays. 
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received  power  (dBm)  normalized  power  (dB) 


(a )  Normalized  power  values 


(b)  Received  power  values 

Figure  20:  Azimuth  profile  for  scene  6  with  two  targets  located  at  -0.4°  and  + 1.6 °, 
processed  for  both  MIMO-1  and  M1MO-2  arrays. 
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received  power  (dBm)  normalized  power  (dB) 


(a )  Normalized  power  values 


(b)  Received  power  values 

Figure  21:  Azimuth  profile  for  scene  4  with  two  targets  located  at  -0.4°  and  +4.6°, 
processed  for  both  MIMO-1  and  M1MO-2  arrays. 
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received  power  (dBm)  normalized  power  (dB) 


(a )  Normalized  power  values 


(b)  Received  power  values 

Figure  22:  Azimuth  profile  for  scene  8  with  two  targets  located  at  + 1.6 °  and  + 14.8 °, 
processed  for  both  MIMO-1  and  M1MO-2  arrays. 
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5  Discussion  and  conclusions 


In  this  report,  we  have  demonstrated  experimentally  that  the  MIMO-1  array  fac¬ 
tor  is  equal  to  the  2-way  phased  array  factor,  as  predicted  with  the  mathematical 
formulation  of  Section  2,  for  uniform  phase  and  magnitude  illumination. 

We  have  further  studied  the  MIMO-2  array  configuration,  for  which  we  used  two 
transmitters  and  seven  receivers.  We  have  demonstrated  that  the  MIMO-2  configu¬ 
ration  provides  improved  resolution  and  much  reduced  complexity,  at  the  expense  of 
gain  loss  and  higher  sidelobes.  The  reduced  gain  can  be  offset  by  the  MIMO  advan¬ 
tage  of  providing  greater  coverage  and  longer  dwell  times  by  illuminating  the  whole 
scene  at  once  with  orthogonal  waveforms.  This  discussion  is  beyond  the  scope  of  this 
report. 

The  higher  sidelobes  of  MIMO-2  array  can  be  circumvented  by  using  appropriate 
weighting  of  the  MIMO  channels  before  summation.  Indeed,  the  MIMO-1  config¬ 
uration  can  be  viewed  as  a  triangle  weighting  of  the  MIMO-2  configuration.  The 
convolution  of  the  transmit  and  receive  arrays,  i.e.  the  virtual  array,  is  the  convolu¬ 
tion  of  two  rectangle  functions  for  the  MIMO-1  configuration  and  results  in  a  triangle 
with  coefficients  [1, 2,  3, 4,  5, 6,  7,  8,  7,  6, 5, 4,  3,  2,  l]/8  in  our  case.  The  virtual  array  of 
the  MIMO-2  configuration  is  an  array  of  14  elements  each  with  a  coefficient  of  1.  In 
Figure  23,  we  have  reproduced  the  results  of  Figure  19  for  scene  2  for  the  MIMO- 
1  and  MIMO-2  arrays.  We  also  have  plotted  the  results  for  MIMO-2  array  using 
a  triangle  weighting  with  coefficients  [1,  2,  3, 4, 5,  6,  7,  7, 6,  5, 4,  3,  2,  l]/8  to  mimic  as 
closely  as  possible  the  MIMO-1  array.  The  results  show  that  using  this  weighting  we 
can  recover  the  results  of  MIMO-1  in  terms  of  beamwidth  and  sidelobes,  with  much 
reduced  complexity,  at  the  expense  of  further  gain  loss. 

As  discussed  in  Annex  A,  an  important  limitation  of  this  experiment  is  the  presence 
of  ground  multipath.  It  has  limited  our  capability  to  demonstrate  the  benefits  of 
increased  resolution  with  the  MIMO-2  array,  given  that  the  two  targets  in  a  scene 
appeared  with  different  power  levels  as  a  result  of  multipath  with  different  heights 
above  the  ground.  In  a  future  experiment,  this  effect  could  be  minimized  by  a  careful 
choice  of  heights,  and  by  performing  the  experiments  over  grass  rather  than  snow  to 
minimize  specular  reflection. 

Another  limitation  of  this  experiment  is  the  time  requirement  to  collect  one  data 
set.  Given  that  we  used  time  division,  one  data  set  required  switching  the  transmit- 
receive  channels  64  times.  A  future  prototype  could  make  use  of  waveforms  with 
orthogonality  based  on  different  bands  or  coded  waveforms,  which  would  require  an 
8-channel  receiver. 
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received  power  (dBm)  normalized  power  (dB) 


(a )  Normalized  power  values 


(b)  Received  power  values 

Figure  23:  Azimuth  profiles  for  scene  2  with  one  target  located  at  +4.6°,  lower 
height,  processed  for  MIMO-1,  MIMO-2  with  and  without  triangle  weighting. 
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Annex  A  Ground  multipath 


In  this  annex,  we  demonstrate  using  a  flat-earth  model  that  ground  multipath  is  a 
likely  cause  of  the  power  level  fluctuations  observed  in  the  experimental  MIMO  data 
presented  in  this  report. 

As  discussed  in  Section  3.1,  the  MIMO  experiments  were  conducted  outside  on  a 
relatively  flat  terrain  covered  by  a  thin  layer  of  snow.  The  terrain  was  flat  but  could 
have  a  small  slope.  This  terrain  can  be  considered  a  good  specular  reflector.  The 
height  of  the  targets  is  known,  relative  to  an  arbitrary  reference  height,  but  the 
absolute  height  relative  to  the  local  ground  is  not  known  precisely.  The  radar  height 
is  estimated  to  be  between  2.3  and  2.45  m  relative  to  the  average  terrain  height. 
Vertical  polarization  was  used. 

We  use  a  flat-earth  model  as  depicted  in  Figure  A.l.  The  relative  height  of  the 
different  targets,  as  measured  with  the  Leica  TS15  Total  Station,  are  tabulated  in 
Table  A.l.  The  targets  were  located  at  a  range  rg  of  45  m  from  the  radar.  As  a  result, 
the  grazing  angle  T  at  the  point  of  reflection  on  the  ground  was  approximately  6°. 
Both  reflected  and  direct  path  rays  are  within  the  main  beam  of  the  antenna,  thus 
there  is  no  antenna  pattern  effect  in  elevation. 


Ta 


Figure  A.l:  Flat-earth  model  for  ground  multipath. 


The  delay  5  between  the  reflected  and  direct  paths  can  be  estimated  as 


Hrefl  If  dir 

C 

=  \  (V {hr  +  htf  +  t\  -  sj ( hr  -  htf  +  r^j 

<  Ins, 


(A.l) 
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Table  A.l:  Targets  relative  height  as  measured  by  the  Leica  TS15  Total  Station. 


Target 

ht 

(m) 

ht  —  hr 
(m) 

PI 

99.9910 

-0.8680 

P2 

100.0233 

-0.8380 

P3 

100.0833 

-0.7780 

P3,  lower  height 

99.5307 

-1.3280 

P4 

100.1721 

-0.6880 

radar 

100.8580 

0 

which  is  much  less  than  the  100  /is  width  of  the  radar  pulse.  We  can  thus  assume 
that  the  reflected  and  direct  path  pulses  interfere  at  the  location  of  the  target.  The 
delay  S  is  also  less  than  the  3  ns  resolution  of  the  radar  (A t  =  c/2B  where  B  is 
the  150  MHz  bandwidth  of  the  radar),  hence  the  different  multipath  signals  are  not 
resolved  by  the  radar. 


Assuming  that  hr,  ht  R,  and  Q\  ~  62  =  0 ,  we  can  use  the  following  formula  for  the 
propagation  factor  F[ 22]: 


f  =  m 


1  +  p2  —  2p  cos 


(A.2) 


where  T  =  is  the  complex  reflection  coefficient  of  amplitude  p  and  phase  0,  A 

is  the  wavelength  and  f{9)  is  the  antenna  pattern  effect.  f{9)  can  be  set  equal  to  1 
in  our  case. 


The  complex  reflection  coefficient  for  vertical  polarization  can  be  determined  using 
the  following  formula  [23]: 

e  sin  T  —  Je  —  (cos  T)2 

Tv  = - \  (A- 3) 

e  sin  T  +  y  e  —  (cos  T) 

where  e  =  e'  —  ie"  is  the  complex  dielectric  constant.  Reported  values  for  the  complex 
dielectric  constant  of  snow  at  9  GHz  range  between  1.15  and  2.6  for  e'  and  between 
0.1  and  0.95  for  e"  depending  on  the  water  content  and  density  of  snow  [24],  The 
reflection  coefficient  around  6°  grazing  angle  converges  to  0.64  for  the  magnitude  and 
178°  for  the  phase  for  most  of  the  reported  dielectric  constant  values.  We  use  these 
values  in  what  follows. 
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In  Figure  A.  2  is  shown  the  propagation  factor  F 2  on  a  logarithmic  scale  as  a  function 
of  the  target  height  relative  to  the  radar  height,  for  an  assumed  radar  height  of 
2.35  m.  The  propagation  factor  is  plotted  for  both  a  smooth  surface  and  surface  with 
roughness  factor  of  0.7.  We  see  that  the  propagation  factor  can  take  values  between 
+4.3  dB  for  constructive  interference  and  -8.8  dB  for  destructive  interference  for  the 
smooth  surface,  or  between  +3.2  dB  and  -5.1  dB  for  the  rough  surface.  Note  that 
the  propagation  factor  F2  appears  in  the  one-way  radar  range  equation,  and  that  FA 
appears  in  the  two-way  radar  range  equation.  The  propagation  factor  can  thus  easily 
explain  some  7  dB  difference  between  the  same  target  at  two  different  heights  (scenes 
2  and  3  of  our  experimental  data). 


Figure  A. 2:  Propagation  factor  as  a  function  of  target  height  relative  to  radar 
height. 


In  fact,  the  dashed  vertical  lines  of  Figure  A. 2  indicate  the  heights  of  the  different 
targets  and  the  radar.  Note  that  we  do  not  know  exactly  the  height  of  the  radar 
relative  to  the  height  of  the  reflection  point  on  the  ground,  so  the  dashed  lines  on 
the  Figure  could  be  shifted  relative  to  the  propagation  factor  curve.  Nonetheless, 
using  our  best  guess  of  radar  height  of  2.35  m  and  the  snow  dielectric  constant,  we 
can  tabulate  the  predicted  values  for  the  propagation  factor  in  the  two-way  radar 
equation.  Table  A. 2  compares  the  predicted  values  with  the  relative  power  of  the 
targets  from  Figure  13.  There  is  a  fairly  good  agreement  between  the  simple  model 
and  the  experimental  results,  except  for  the  target  at  15°  azimuth.  But  recall  that 
this  target  is  losing  an  extra  6  dB  due  to  the  azimuth  antenna  beamwidth.  The 
predictions  with  a  smooth  surface  exacerbate  the  difference  in  power  levels,  where 
the  predictions  using  a  rough  surface  better  compare  with  the  experimental  results. 
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Table  A. 2:  Power  levels  predicted  by  the  simple  flat  earth  model  and  observed  in 
the  experimental  data. 


Observed 

Predicted  two-way 

relative 

propagation  factor 

power 

smooth  surface 

rough  surface 

Target 

(dB) 

(dB) 

(dB) 

PI 

44.9 

7.1 

5.1 

P2 

43.3 

4.0 

2.5 

P3 

33.8 

-10.3 

-7.4 

P3,  lower  height 

40.0 

-0.2 

-1.0 

P4 

33.0 

-0.1 

-0.9 
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